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0. Introduction and Preliminaries 



Electronic spectroscopy of partly-filled shell ions, like and ions, embedded in con- 
densed matter surroundings has been the object of important developments in the recent 
years. In particular, the advent of tunable dye lasers has made possible to perform two- 
photon spectroscopy experiments. This Letter is concerned with p-photon spectroscopy of 
a transition ion with ground atomic configuration n£^ in a given crystal field with generic 
symmetry G. We address here the question of obtaining a closed-form expression for the in- 
tensity of a p-photon (absorption) transition between an initial Stark level i with symmetry 
species F and a final Stark level / with symmetry species F' . The labels F and F' stand for 
two IRC's (irreducible representation classes) of the double group G* of G ~ G* jZ^. The 
group G* may be considered (up to isomorphism) as a subgroup of SU{2). 

According to Dirac and Goppert-Mayer, the transition moment Mj(-r^^_,j(r'7') ^or a 
p-photon transition between the initial state i{T^) and the final state f{T'Y) is given by 

M,(r^)_^(r.y) = Yl (^^^^^1 l^i) I 1^^) - • - (^^-1 1 1^^^) + perm (1) 

Vj 

in the framework of the electric dipolar approximation. The sum on Vj in (1) has to be 
performed on all the intermediate states and "perm" indicates that p\ permuted terms arise 
in the Feynman representation associated to Mj(r7)^/(r'7') when the p photons are different. 
The operator (D. E/.) in (1) is the scalar product of the dipole moment operator D for the N 
electrons with the unit polarization vector of the single-mode for the kth photon. (We use 
single-modes to describe the electromagnetic field.) Furthermore, the factors (flj — ily.) are 
energy denominators for the system radiation field plus ion in its environment. Finally, the 
labels 7 and 7' distinguish the components of the Stark levels i and / when the dimensions 
[F] and [F'] of the IRC's F and F' are greater than 1, respectively. 
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At this stage, we must consider two types of p-photon transitions: the parity- allowed 
and parity-forbidden transitions. The parity-allowed transitions occur either between levels 
of the same parity {i e n£^ and / e n£^) when p is even or between levels of opposite 
parities (i e n£^ and / e n£^~^n'i' with (—1)^+^' = —1) when p is odd. The parity- 
forbidden transitions {n£^ n£^ for p odd and n£^ n£^~^n'£' with (—1)^+^ = —1 for 
p even) may become possible, via the Van Vleck mechanism, if G does not have a center of 
inversion (or if the inversion symmetry is broken by vibronic coupling or ligand polarization) . 

The sums over the intermediate states vj in (1) can be effectuated by using a quasi- 
closure approximation of the type of the ones used by many authors [1-5] for p = 1 or 2 in 
the cases of and ions. Such an approximation used in conjunction with recoupling 
techniques for SU (2) irreducible tensorial sets leads to 

Mi(r^)^^(r,y) = {fT'y\H,s\iT^), (2) 

where the model-dependent operator H^fi is an effective operator described in Section 2. 

Both for parity-allowed and parity-forbidden transitions, we take the initial and final 
state vectors in the symmetry adapted form (a; = z or /) 

|a;r7) = J2 laSLJaT^) c{aSLJaT;x). (3) 

aSLJa 

The initial state i{T^) belongs to the configuration n£^ and the final state /(r'7') belongs 
either to the configuration n£^ or to the configuration n£^~^n'£' . In (3), the SU{2) D G* 

vectors of type \ jaT^) are connected to the SU{2) D U{1) vectors of type \ jm) by 

j 

\jar^) := ^ \jm) {jm\jar^). (4) 

m=-j 

The branching label a in (4) has to be used when the IRC F of G* occurs several times in 
the IRC (j) of SU{2). This external multiplicity label may be described, at least partially, 
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by IRC's of subgroups of SU{2) that contain in turn the group G*. The expansions (3) and 
(4) are of a completely different nature. Equation (4) is of a group theoretical nature: the 
reduction coefficients {jm\jaT^) depend on the chain SU (2) D G* but not on the physics 
of the problem. On the other hand, the expansions (3) depend on physical models: the c 
coefficients may be obtained from the diagonalization of one or two G-invariant Hamiltoni- 
ans (according to whether as the initial and final configurations are the same or different) 
describing the ion in its environment. It is to be noted that, thanks to Schur's lemma, we 
may always standardize the reduction coefficients {jm\jaTj) in (4) in such a way that the 
expansion coefficients c in (3) be independent of the generalized magnetic quantum numbers 
7 and 7'. This result is essential in view of the sum on 7 and 7' in (5) below. Equation (3) 
indicates that we adopt here the philosophy of [6-9] where the diagonalization-optimization 
of the G-invariant Hamiltonian(s) is done, without loss of generality, in a weak-field basis 
adapted to the chain SU{2) D G*: the labels a, S, L and J are not taken a priori as good 
quantum numbers; the only good quantum numbers are r,7 and r',7' for the initial and 
final state vectors, respectively (and the parity for parity-allowed transitions). 

The physical quantity we want to calculate in Section 2 is the intensity 

Si{r)^f{r') ■■= Yl l^i(r7)-/(r'7')P (^) 

77' 

for the i(r) — > /(F') transition. The sum in (5) has to be done over all the external 
multiplicity labels 7 and 7', i.e., over all the components of the initial and final Stark levels. 
(The labels 7 and 7' may be characterized, at least partially, by IRC's of subgroups of G*.) 
The purpose of this work is to show how far group theory can go in the determination 
of 5'j(-r)^/(r') in connection with sophisticated models for the operator i^eff of (2). In 
particular, we want to clearly separate the elements arising from group theory (via Wigner- 
Racah calculus) from the model-dependent quantities (both for the light-matter interaction 
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and the description of the £^ ion in its environment). Rather than considering G* as an 
isolated group, we shall consider it as a subgroup of SU{2), an approach that has already 
proved to be very useful as far as the energy spectrum of the £^ ion in G symmetry is 
concerned [8, 9]. The basic tool is thus the Wigner-Racah algebra of the group G* in a 
G* C SU{2) basis. Therefore we develop, in Section 1, those nontrivial aspects of the 
Wigner-Racah algebra for the chain SU{2) D G* that are essential for performing the sum 
in (5). The corollary obtained in Section 1 is applied in Section 2 to derive the sum rule 
(23) with the accompanying properties and rules (26-28). 



1. Sum Rule for SU{2) D G* Coupling Coefficients 



Equation (4) defines the {jaV'^} scheme for the chain SU{2) D G* that is more appropriate in 
condensed matter spectroscopy than the {jm} scheme for the chain SU{2) D U{1) [6, 8, 9]. 
For j, a and F fixed, the set {\jaT^) : 7 = 1, • • • , [F]} is a (standardized) G-irreducible 
tensorial set (in the sense of Fano and Racah) of vectors associated to F. Similarly, from 
the SU{2) D U{1) spherical tensor operators Uq^\ we define the operators by an 

expansion similar to (4) so that, for /c, a and F fixed, the set |?7^r7 • 7 = ' " " ; | is a 
G-irreducible tensorial set of operators associated to F. The latter G-irreducible tensorial 
sets are also labelled by IRC's of SU{2) and, therefore, we can easily generate, by direct 
sum, nonstandard SU (2)-irreducible tensorial sets. Thus, we may apply the Wigner-Eckart 
theorem for the group SU{2) in a nonstandard basis adapted to its subgroup G* . As a 
result, we have [8] 

(rijiaiFi7i|t/^P^|T2j202r272) = {riii\\U^^^\\T2j2) 

El I /( " ' " ^' 

03r373 V a3r373 aiFi7i / V a3F373 aF7 a2F272 




(7) 



with 

/ ji h h \ 1 3i h h 

\airi7i a2r272 a^T^^^ J mimams \mi 1712 ms 

UimiljiaiTi^i)* (j2m2|j2a2r272)* {jsmsljsasT 3^3)' 

and 

j \ f j j \ 

- [J]'^' 7 |, (8) 

airi7i a2r272/ \a1r171 aoro7o a2r272 

where Fq denotes the identity IRC of G*. In the right hand side of (6), the quantum numbers 
Ti and T2, external to the chain SU{2) D G*, appear in the reduced matrix element. 

The / or S-jaT^ symbol in (7) is an SU (2) D G* symmetry adapted form of the usual 
S-jm Wigner symbol. It constitutes a symmetrized form of the coefficient [8] 



(9) 



/ h j2 k \ ( ji k j2 

\airi7i a2r272 ClTI J miQm2 \— TOi q 777.2 

{jimi\jiaiTi-fiy (kqlkaT'j) {j2m2\j2a2^2l2) 
that generalizes the / coefficient defined by Racah and some of his students (in particular 

Schonfeld and Flato, see references [6, 7]). 

The 2-jaT^ symbol defined by (8) turns out to be an SU (2) D G* symmetry adapted 
form of the Herring- Wigner metric tensor. Indeed, from (7) and (8) we have 
j 

airi7i a2r272 

Therefore, the metric tensor given by (10) allows us to handle all the phases occurring in 
the {jaV^} scheme [8]. 

By combining (7), (9) and (10), we can rewrite (6) in the simple form 



= J2 (-1)'+"^ (jm|jairi7i)* (j,-m|ja2r272)*. (10) 



(rijiairi7i|t/i? Jr2i2a2r272) = (riji||t/('=)||r2j2) / • (11) 



ji 32 k 

airi7i a2r272 ar7 



The interest of (6) and (11) for electronic spectroscopy of ions in crystalline fields has been 
discussed in [8, 9]. From a mathematical viewpoint, it is to be observed that the factorization 
in (11) into the product of a G*-dependent factor (the / coefficient) by a G*-independent 
factor is valid whatsoever the group G* is multiplicity free or not. The internal multiplicity 
problem arising when G* is not multiplicity free is thus easily solved by making use of (11). 

The / coefficients may be considered as the components of a third-rank tensor (second- 
rank covariant and first-rank contravariant). In contradistinction, the / symbol defines a 
third-rank contravariant tensor. The other properties of the / and / coefficients that are of 
relevance for this work are based on the following lemma. 

LEMMA (Racah's lemma). We have the factorization formula [10] 



where {j20-2^2-\-kaV\jiaif3Vi) is an isoscalar factor or reduced f symbol (independent of the 
quantum numbers 71, 72 and ^) and (r2r727|r2r/3ri7i) a Clebsch-Gordan coefficient for 
the group G* considered as an isolated entity. (The label j3 is an internal multiplicity label 
to be used when the Kronecker product T2 <S>T is not multiplicity free.) 

Proof. Equation (12) is a consequence of Schur's lemma. 

In order to be able to use tables and/or computer programs for the highly symmetrical 
/ symbol (and for the i?/3 or reduced f symbol, see (13) and (14)), it is useful to transcribe 




(12) 



J2 iha2r2 + kaT\jia^pT^y (r2r727|r2r/?ri7i)*. 
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(12) for the / symbol. We thus obtain 

\oiri7i a2r272 a3r373 / /3 \airi a2r2 asFs / \7i 72 73 

(13) 

where Vg is an extension, taking into account the multiplicity label /3, of the V coefficient 
introduced by Griffith [11]. The coefficient 

31 32 3z \ / Ji 32 33 \ 

Rpl ]=Rp[ (14) 

\airi a2r2 asTs J VaiFi a2r2 ) q* 

characterizes the restriction from SU{2) to G*; it is a generalization of the V symbol intro- 
duced by the Chinese school of Tang Au-chin [12]. (In (13) and (14), the label /3 is absolutely 
necessary when the identity IRC Fq appears several times in Fi (8)F2 (8)F3.) The Rp symbol 
is connected to the isoscalar factor of (12) by 

(i2a2r2 + /carijiai/^iFi) = (-1)^^ [T^]-' ^ J] UiT.T^T)^^,,^, 

k 32 ii \ / ii ii \ * (15) 



aV a2F2 a^Ts / \aoTo aiFi a3F3 

where U (F1F2F) is an arbitrary unitary matrix. Indeed, this matrix is a simple phase factor 
when the Kronecker product F2 (8) F is multiplicity free. 

Repeated applications of Racah's lemma and of the orthonormality-completeness prop- 
erty for the Clebsch-Gordan coefficients of G* lead to the following corollary. 

COROLLARY. The sum 

( 31 32 k3 \* ( j4 3b h \ 

^-E/ / (16) 

7172 V aiFi7i a2F272 a3F373 / \ aj"i^i a5F272 a6r676 / 



is diagonal in {V^^Vq) and (73,76)- Furthermore, the diagonal value of X is independent of 
73. More precisely, we have 

X = 5{T,,T,) 5(76,73) (-l)2(^3+fca) (f^-^] [r^]-! [r^] 

(17) 

^(iioi, ^2^2, ^303,^404,^505, A;606; rir2r3), 

where 

■^(jioi, ^202, A;3a3,j4a4, ^505, Ajeae; rir2r3) := 

(j202r2 + A;303r3|jiai/3ri) (j505r2 + fc606r3|i404/3ri)* ^^^^ 

13 

that reduces to the product of two isoscalar factors when G* is multiplicity free. 

It is also useful to express Z and thus X in terms of the symbols which are more 
symmetrical than the isoscalars factors. In this respect, by using (15) we can write (18) as 

^(jlOi,j202,A;303,i4a4,i505,A:606;rir2r3) = (-l)2('=3+fc6) [r^]-2 ^-^j 

/3 a^r^ ayVy \aoTo aiTi a^Tj: / \a3T3 a2T2 a^T^ / (19) 

(0 j4 j4 \ / ke j5 j4 

aoTo a4Ti ayVy J \ aoTs a5T2 ayVy 

that is in fact independent of the U matrix of Equation (15). 



2. Intensity Formula 

For a given transition operator H^fi, we are now in a position to effectuate the summation 
over 7 and 7' in (5) and, thus, to obtain a compact expression for the intensity 5'j(r)^/(r')- 
A very general expression for the effective operator i^eff is 

^eff= E EE E^[^i^2---Vi;^;^s/cl;/c] (^{p('=p-i)xW}«.w('=^^-)'=), 

fcife2-"fep_i t kskL k 

(20) 
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where ( . ) is a scalar product involving electronic (W^'^'^'^^^'^) and nonelectronic (p('^p-i) 
and X*^*-*) tensors. In Equation (20), p('^p-i) denotes the polarization tensor defined by 

p(fcp-i) :={... {{EiE2}^''^)E3}('=2) . . . Ep}('=''-i) (21) 

that describes the coupling of the p polarization vectors and that is entirely under the 
control of the experimentalist. The tensors 'W('^skL)k ^^^^ ^(t) relative to the ion and 
its environment: W^'^'^'^^^'^ is an electronic double tensor of spin degree ks, orbital degree 
kL and total degree k whereas X^*) is a (single) tensor of the degree t, for instance, the 
ligand polarization tensor or the crystal-field tensor. (The tensor X^*) may result of the 
coupling of several G-invariant operators according to the order of the mechanism used.) 
The C coefficients are model-dependent parameters that depend, among other things, on 
the initial, final and intermediate configurations as well as on the energies of the p photons. 
Equation (20) unifies various models introduced in [1-4] for p = 1 and in [5, 13-24] for p = 2. 
It can be applied to both (i) parity- allowed and (ii) parity-forbidden transitions. 

(i) The particular case of intra-configurational two-photon transitions {p — 2) is of 
special interest. For parity-allowed transitions within the configuration n£^ , it is sufficient 
to take t = with X(°) = 1 in order to describe second-order mechanisms when ks — 
and second- plus third-order mechanisms when ks = and ks ^ 0- The situation where 
t ^ permits to describe (other) third-order mechanisms when ks = and fourth-order 
mechanisms when A;^ 7^ 0. As a matter of fact, (20) gives back the model introduced by 
Axe [5] (with t = 0, ^5 = 0) for second-order mechanisms and extended by different authors 
[14, 15, 17, 18, 20, 21] (with t 0, ks ^ 0) to take into account higher-order mechanisms. 

(ii) Two particular cases deserve a special attention, viz., the intra-configurational 
one-photon transitions and the inter-configurational two-photon transitions. For parity- 
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forbidden transitions, X^*) is either the hgand polarization tensor or the odd crystal-field 
tensor. In the case of one- photon transitions {p = 1) within the configuration n£^ , (20) 
with ks = gives as a special case the operator implicitly considered by several people 
in order to explain intra-configurational one-photon transitions for [1, 2] or [3, 4] 
ions in crystals. In the case of inter-configurational two-photon transitions {p = 2) between 
the configuration n£^ and the configuration n£^~^n'£' with £ + £' odd, (20) with ks = 
yields as particular cases the operators introduced in connection with either a static cou- 
pling mechanism [13, 16, 19, 22, 24] when X^*) is the odd crystal-field tensor or a dynamic 
coupling mechanism [19, 22] when X*^*) is the ligand polarization tensor. 

The transition matrix element for the operator (20) may be easily calculated owing to 
(11). We thus get 



M,(r^)^^(r'y) = E E (-^V'"' c{a' S' L' J' a'V ; /)* c{aSLJaT; i) 

a'S'L'J'a' aSLJa 

-it) 

370 

kik2---kp-i t kskL k ao 

(22) 



E E E E E ^ [^1^2 ■ ■ ■ /cp-i; t- kskL-, k] x^^^, 

;2'"fep_i t kskL k ao 



k kp-i t \* ( J J' k 

f\ / 

h"V"i' a"V"i' aoTo^oJ \aT^ a'T'Y b"T"Y' 



(t) 

do To 70 



where we have used the fact that the tensor X*^*) has only G-invariant components X 
distinguishable by the multiplicity label ao- The right hand side of (22) can be expressed 
in terms of / coefficients in view of the connecting formula (arising from the comparison of 
(6) and (11)) between the / and / symbols. 

By introducing (22), in terms of / or / symbols, into (5) and by using the corollary of 
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Section 1, we obtain the intensity formula 

^^(r)^/(r') = E E E E E ^[{^^} rsT"; TV'] p!:'^'' (^^-y'T (^3) 

r s r" 7" 

(with 1 < i < p — 1), the form of which holds for both parity- allowed and parity-forbidden 
p-photon transitions. The intensity parameters / in (23) are given by 

/P.}{«r.r";rn := EEEEEE 

J 'a' J a J 'a' J a kr is 

YkiJ'a'V, JaV, {ki} rfV") Ye{J'a'r', JaV, {£i} ssT")* (24) 

Z{Ja, J' a, kf, Ja, J' a', £s; FFT"), 
where the Y symbol is defined through 

YkiJ'a'r',jar,{k,}rfr") mnv'^" nm)'^" (-i)'-"' E E E 

a'S'L' aSL ksk^ 

c{a'S'L'J'a'r';f)*c{aSLJar;i) (a'S'L'J'\\W^''^''^^''\\aSLJ^ ^^5) 
J2 C [k,k2 ■ ■ ■ k^^,;t; kskL-, k] Xi%^^ f( 

tao \fT"y' rV"i' aoroTo. 

which does not depend on 7". The / parameters in (24) can be calculated in an ah initio way 
or can be considered as phenomenological parameters. In both approaches, the following 
properties and rules are of central importance. 

PROPERTY 1. In the general case, we have the (hermiticity) property 

I[{ei}{h}srT";TT'] = /[{/cj rsF"; (26) 
that ensures that /S'j(r)_»/(r') is a real number. 

PROPERTY 2. In the case where the group G is multiplicity free, we have the factorization 
formula 

I[{h}{QrsV"-VV'] = x[{h}rV"-VV'] xmsV"-VV']\ (27a) 
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where the function x is defined through 

x[{ki}rV"-VV'] =^Y.Y. Yk{J'a'V',JaV,{ki}rrV") {J' a'V + kfV"\JaV) (276) 

J'a' Ja kf 

which follows from (24) and the factorized form of Z (see (18)). 

The number of independent parameters / in the expansion (23) can be a priori deter- 
mined from the two foUowing selection rules used together with Properties 1 and 2. 

RULE 1. In order to have 'S'j(r)-«-/(r') 0? necessary that 

T" c r'* r r" c {kp-i) v" c (£p_i) r" c {k) v" c (^), (28) 

where (kp-i), {£p-i), (k) and {£) are IRC's of the group 0(3) and T'* is the complex conju- 
gate of v. 

RULE 2. The sum over {ki} and {£i} in the intensity formula (23) is partially controlled 
by the following points, (i) The order of the mechanism used for describing the absorption 
processus: kp-i and £p-i cannot vanish if only p- order mechanisms, corresponding to t = 
ks — 0, are taken into consideration. (In this case, kp^i = k^ = k and £p-i = £l = £■) 
Conversely, kp-i and £p-i may be zero if higher- order mechanisms, corresponding to t ^ 
and/or ks ^ 0, are introduced, (ii) The nature of the photons: if the p absorbed photons 
have the same polarization, then ki{and £i) = 0,2 so that if ki{or £i) = we have ki+i{or 
£i+i) = 1 and A;j+2(or £i+2) = 0, 2 when i > 1. 

3. Discussion and Closing Remarks 

For low values of kp-i and £p-i, there is no summation on r and s, two branching multiplicity 
labels of type a, in the intensity formula (23). (The frequency of T" in (kp^i) and {£p-i) 
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is rarely greater than 1 for p < 2.) The group-theoretical selection rules (28) impose strong 
limitations on the summation over F" in (23) once F and F' are fixed and the range of values 
of k and £ is chosen. 

The number of independent intensity parameters / in (23) is determined by: (i) the 
number of photons and their nature (polarization, energy), cf. Rules 1 and 2; (ii) the symme- 
try group G, cf. Rule 1; (iii) the symmetry property (26), cf. Property 1; (iv) the use of t = 
and ks = (p-order mechanism) or t 7^ and /cg 7^ (higher-order mechanisms), cf. Rule 2; 
(v) the kind of the (weak-, intermediate- or strong- field) coupling used for the state vectors, 
cf. (3). Points (i)-(ii) depend on external physical conditions. On the other hand, points 

(iv) and (v) are model-dependent. In particular, in the case where the J-mixing, cf. point 

(v) , can be neglected, a situation of interest for lanthanide ions, the summations on k and £ 
in (24) are reduced by the rule \ J — J'\ < k,£ < J + J' , where J and J' are the total angular 
momenta for the initial and final states, respectively. Similar restrictions apply to ks and 
kL in (25) if the S- and L-mixing are neglected. 

The computation, via Equations (18) or (19), (24) and (25), of the / parameters is 
a difficult task in general. Therefore, they may be considered, at least in a first step, 
as phenomenological parameters. In this respect. Equations (18) or (19), (24) and (25) 
should serve as a guide for reducing the number of / parameters. Once the number of 
independent parameters / in the intensity formula (23) has been determined, we can obtain 
the polarization dependence of the intensity strength S'j(r)->-/(r') by calculating the tensor 
products -Pa"r"7" (with K = kp-i,£p-i and a" = r, s) occurring in Equation (23). 

To close this Letter, some remarks are in order. In the particular case p = 1 and 2, the 
Hamiltonian model given by (20) and (21) unifies various models described in the literature 
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for rare earth ions and transition-metal ions. The originahty of this work rests on the use of 
symmetry adaptation methods for the chain SU (2) D G* in conjunction with a very general 
Hamiltonian model for describing simultaneous absorption of p photons between Stark levels 
(rather than between J levels). Nontrivial aspects of symmetry adaptation methods have 
been taken into account in a quantitative way: factorization d la Wigner-Eckart for the chain 
SU{2) D G*, Racah's lemma for the corresponding coupling coeflBcients and orthogonality- 
completeness of the latter. As a net result, we have obtained the intensity formula (23) 
where the dynamics appears in the / parameters and the geometry is contained in the P 
factors describing the polarization dependence. The intensity formula (23) for multi-photon 
absorption can be extended to multi-photon emission as well as to Rayleigh and Raman 
scattering modulo some caution with the C parameters in (25). The general form of (23) is 
also valid for other multi-photon processes, as for example the simultaneous absorption of 
several photons, certain by electric-dipole absorption and others by magnetic-dipole and/or 
electric-quadrupole absorption. 

A word should be said about previous works on this subject. The subject treated in 
the present Letter has been touched upon in [13, 21, 23-26] for p = 2 and fully developed 
in the thesis by one of the authors (M.D.) for p arbitrary. The accent has been put here 
on the intensity, rather than on the transition moment as in some works dealing with two- 
photon spectroscopy [13, 21]. Furthermore, this work is concerned with models and their use 
in connection with symmetry adaptation methods, rather than with qualitative symmetry 
considerations only. In this respect, our work represents a further important step besides 
the pioneer works by Inoue and Toyozawa [25] , on one hand, and by Bader and Gold [26] , on 
the other hand, where only symmetry considerations, arising from the group G considered 
as an isolated entity, are introduced in the situation where p — 2. 
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